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Abstract 

The model under consideration is a semi-infinite two-dimensional two-component 
plasma (Coulomb gas), stable against bulk collapse for the dimensionless coupling 
constant /3 < 2, in contact with a dielectric wall of dielectric constant = 0. The 
model is mapped onto an integrable sine-Gordon theory with a "free" Neumann 
boundary condition. Using recent results on a reflection relationship between the 
boundary Liouville and sine-Gordon theories, an explicit expression is derived for the 
surface tension at a rectilinear dielectric - Coulomb gas interface. This expression 
reproduces the Debye-Hiickel /3 ^ limit and the exact result at the bulk collapse 
border, the free-fermion point /3 = 2, where the surface tension keeps a finite value. 
The surface collapse, identified with the divergence of the surface tension, occurs at 
P = 3. 

KEY WORDS: Two-component plasma; two dimensions; boundary sine-Gordon model; 
surface tension. 
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1 Introduction 



In a previous paper IQ, the bulk thermodynamic properties (free energy, specific heat, 
etc. . . ) of the two-dimensional (2D) two-component plasma (with logarithmic interac- 
tions), or Coulomb gas, have been obtained exactly in the whole temperature range for 
which the point-particle model is stable against collapse, i.e., for the dimensionless cou- 
pling constant (inverse temperature) (3 <2. The mapping onto a bulk sine-Gordon field 
theory was made and recent results about that field theory 0, @, were used. In a next 
paper a surface property of the same model in contact with an ideal conductor of 
dielectric constant oo, impenetrable to particles, was considered. In particular, the 
surface tension at a rectilinear conductor - Coulomb gas interface was obtained exactly 
as a function of the bulk density, the applied potential and the temperature. A mapping 
onto an integrable boundary sine-Gordon theory [^] with a Dirichlet boundary condition 
was established, and known results about that field theory 0,0, were used. 

The model considered in the present paper is the 2D Coulomb gas in contact with an 
ideal dielectric wall of dielectric constant = 0. The system is mapped onto an integrable 
sine-Gordon theory, now with a "free" Neumann boundary condition. Using very recent 
results on a "reflection" relationship between the boundary Liouville and sine-Gordon 



theories |g[ - ||TOl , an explicit expression is derived for the surface tension at a rectilinear 
dielectric - Coulomb gas interface as a function of the bulk density, or the fugacity, and the 
temperature. The surface tension is checked on its high-temperature expansion derived 
from a renormalized Mayer expansion and on the exact result at the bulk collapse border. 



the free-fermion point /3 = 2, where it keeps a finite value ||TT|. The divergency of the 
surface tension occurs at /5 = 3. 

The model under consideration mimics the interface between an electrolyte (the 2D 
two-component plasma, made of two species of point particles, of opposite charges ibl) 
and an ideal dielectric wall. Classical equilibrium statistical mechanics is used. In the 
grand-canonical formalism, the control parameters are the inverse temperature (3 and 
the two fugacities 2;+ and of the positive and negative particles, respectively. Due to 



the charge neutrahty, the bulk properties of the plasma depend only on the combination 
z = {z+z^Y^'^ ||l2i- It is a specific feature of the present model that also its surface 



properties depend only on z. 

In the bulk of the plasma, the interaction Boltzmann factor of a positive-negative pair 
of charges, r~^, is integrable at small r if and only if /5 < 2 (at large r, the interaction is 
screened), with (3 = 2 being the bulk collapse border of the plasma. In the region (3 > 2, 
one has to attach to particles a small hard core a in order to prevent the collapse. As has 
been argued long time ago |T^, in the limit a ^ 0, while the free energy and the internal 
energy per particle diverge, the specific heat, truncated correlation functions, etc. remain 
finite. In this sense one should also consider our exact result for the surface tension when 
(3 > 2. The stability range for surface properties of the model is (3 < 3. The surface 
collapse at /3 = 3 is identified with the divergency of the surface tension. 

Series expansions and exact results are rather rare for the present 2D configuration 
electrolyte - dielectric wall. More success was attained in the case when the electrolyte 
is modeled by a one-component plasma, i.e., a system of moving point particles of the 
same charge embedded in a rigid neutralizing background. The (3^0 limit was treated 
in ref. Ijl^ and the (3 = 2 free fermion case was solved by Smith (thermodynamics) 
and Jancovici ||16| (density profile, pair correlations, sum rules). When the electrolyte is 
modeled by a two-component plasma, at (3 = 2, the exact solution was available only for 



a hard wall and an ideal conductor wall [jT^, [T^, which admit a Thirring representation 
with an inhomogeneous mass going to and to oo within the wall, respectively. Only 
recently ||lT|, an exact solution has been obtained also for the dielectric wall of interest 
by the method of pfaffians, and that solution offers via an explicit density profile a check 
of the surface tension formula at (3 = 2. 

One should emphasize that the present result was obtained due to recent progress 
in integrable field theories, in particular, a refiection relationship between the boundary 
Liouville and sine-Gordon theories [Q - |10|. Mainly the work |T^, providing an exact 
formula for the one-point function of an exponential boundary operator in the boundary 
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sine-Gordon model with a general integrable boundary action, was crucial. The sur- 
face tension is obtained as a result of an integration from the conductor wall (Dirichlet 
boundary condition) with a known formula for the surface tension |Q, to the dielectric 
wall (Neumann boundary condition), via a continuous class of integrable boundary field 
theories. 

The paper is organized as follows. In Section 2, the mapping of the Coulomb gas in 
contact with an ideal dielectric onto a sine-Gordon field theory with a Neumann boundary 
condition is made. The reflection relationship between the boundary Liouville and sine- 
Gordon theories is briefly explained and basic formulae are written down in Section 3. 
The exact formula for the surface tension is derived in Section 4. Its high temperature 
expansion and the (3 = 2 solution are checked in Section 5. Concluding remarks about a 
possible mechanism for the surface collapse are given in Section 6. 

2 Mapping 

We consider an infinite 2D space of points r G -R^ defined by Cartesian coordinates (x, y). 
The model interface is localized at x = 0, along the y axis. The half-space x < 0, 
impenetrable to particles, is assumed to be occupied by an ideal dielectric of dielectric 
constant = 0. The electrolyte in the complementary half-space x > is modeled by the 
classical 2D TCP of point particles {j} of charge {qj = ±1}, immersed in a homogeneous 
medium of dielectric constant = 1. Classical equilibrium statistical mechanics is used. 
In the grand-canonical formalism, the control parameters are the inverse temperature /3 
and the fugacities z+ and Z- of the positive and negative particles, respectively. For 
the dielectric wall of interest, only strictly neutral charge configurations survive and, 
therefore, the thermodynamics depends only on z = [zj^zJ]^!'^ as will follow directly from 
the formalism presented here. Due to the translational invariance in the y direction, the 
particle densities n^ir) = n_{r) depend only on x. Let us denote their asymptotical 
X oo values hj = n_ = n/2 where n is the total particle number density. 
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The interaction energy E of particles {qj, Vj = {xj > 0, yj)} consists of two parts (see, 
e.g., H), 

(1) 



i<j 



where r* = {—x,y) and v{r) = — ln(|r|/ro) is the 2D Coulomb potential (ro will be set 
for simplicity to unity). The first term corresponds to direct particle-particle interactions, 
while the second one to interactions of particles with the images due to the presence of 
the dielectric wall. Introducing the microscopic charge + image charge density 



Pi^) = J2 Qj [^{^ -Xj) + 6{x + Xj)]6{y - yj) 
j 

the energy (|1|) can be written as 

E = i y dVy dV p{r)v{r, r')p{r') - ^Nv{0) 



(2) 



(3) 



where the integrations over r and are taken over the whole 2D space and — (A^/2)f (0) 
is the self-energy. 

The grand-canonical partition function is defined by 



Z\ Z 



Q{N+,N. 



(4a) 



where 



N 



Q(iV+, iV_) = I n d'r, exp [-/5E({g„ r,})] (4b) 

is the canonical partition function of positive and negative charges and = 
+ A^_. Since, in infinite space, — A/(27r) is the inverse operator of v{r), the Boltzmann 
factor of the interaction energy is expressible as 

/ V(f) exp [/ dV (i0A0 + 



exp 



4 



d r d r' p{r)v{r, r')p{r'] 



(5) 



/ V(f)exp J d^ri 

where (j){r) is a real scalar field and / denotes the functional integration over this field. 
Inserting p from (0), one recognizes in the action of the field theory (^ a nonlocal term 
i-\/TiPJ2j Qj[(p{^j^ Uj) + 4>{~^jj Uj)]- To make the field theory local, we shall reformulate it 
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as a boundary problem by introducing two new fields 



<Pe{x,y) = -j=[(l){x,y) + (f){-x,y)] 



(6a) 
(6b) 



defined only in the positive x > half-space. The even field has a Neumann boundary 



condition dx(pe\x=o 
Since it holds 



and the odd field has a Dirichlet boundary condition 0o 

dV {(f)eA(j)e + (j)oA(j)o) 



x=0 



dV- 



1 

2 Jx>o 



0. 



(7) 



the odd field, contributing only by its free-field part (j)oA(po/'2, disappears from (||) by the 
numerator-denominator cancelation. By integration per partes, the term 0eA0e can be 
rewritten as — (V0e)^, with a vanishing contribution from the boundary. The rhs of (^) 
is thus expressible as 

/P0eexp{4>odV [-|(V0e)'] +iyWE,g,0e(r,)} 



(8) 



After inserting (|]) via (^ into (^, one proceeds along the standard line and finally 
expresses S in terms of the 2D Euclidean sine-Gordon theory formulated in the half-space 

X > 0: 

f V(j)exv ( L.n dV [-ifV(i))2 + 2z cos^^/2i;50)] } 



/P0exp [-/^^,d2ri(V0)^ 
with Neumann condition dx4>{x,y)\x=o = for the field at the boundary. Here, z is the 
fugacity (z^z^Y^'^ renormalized by the self-energy term exp[/3t'(0)/2], and the uniformly 
shifted 0e is renamed as 0, 

\n{z+/z_) 



+ 



(10) 



The dependence of the statistics exclusively on (2;_|_2;_)^/^ is now evident. The rescaling of 
the field d) 6/ \/2n transforms (0) into the form 



JVcj) exv[-AsG{z)] 



(11a) 
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with the action 

1 



Asg{z) = / dV 

Jx>0 



-^(V0)^ -22cos(260) 
An 



(lib) 



h = \fP (11c) 

and boundary condition dx(t)\x=o = 0, which is more convenient for our purpose. Here- 
inafter, the dependence of quantities on the temperature parameter h will be omitted in 
the notation. 

The boundary sine-Gordon theory (0) is a member of integrable sine-Gordon theories 
defined in the half-space x > [§: 

^^'^'^^ - JV<P exp[-AsG{z = 0,zb)] ^''^^ 

with the action 

/' r 1 1 

Asg{z,zb)= dV — (V(/))2 - 22cos(26</)) - 2zb dy cos{b(f)B) (12b) 

Jx>0 14:71 J J-OO 

Here, 4>B{y) = <P{^ = 0,?/) is the boundary field and zb the boundary fugacity. The 
underlying sine-Gordon theory (0) with Neumann boundary condition dx4>\x=o = is 
known to correspond to the "free" case Zb = (see, e.g., P, ||21|). The model of the 
metal-electrolyte boundary studied in ref. 0] corresponds to the limit zb oo, fixing the 
value of the boundary field 0|x-o = (Dirichlet boundary conditions for zero potential 
difference between the metal and the electrolyte). 

The grand potential Q = —(3^^ lnS(2;, zb) is the sum of a volume part and a surface 
part: 

n = -Vp{z) + S^{z,ZB) (13) 

where p is the bulk pressure, dependent only on z, and 7 the surface tension, dependent 
on both z and zb- For a strip L x /2, i? 00 in the y direction and L large in the x 
direction, the "specific" Vl/Ris given by 

lim ^ = -Lp{z) + 7(2;, Zb) (14) 



Thus, with regard to (jT^) and taking into account the crucial z^-independence of p, one 
finds 



(15) 



where the obvious symmetry relation 



(16) 



implied by the invariance of the action (p.2b|) with respect to the transformation 



was assumed. (. . ■)z,zb denotes the averaging with the sine-Gordon action (|12b|). 



3 Reflection property 



The Liouville field theory, formulated in the half-space x > with a conformally invariant 
condition at a; = 0, is defined by the action 

1 



dV 



x>0 



47r 



(V0)^ + ^e^'"^ 



^Zb\ dy e" 



(17) 



and the boundary condition on the field at infinity, (f){x,y) = — (51n(x^) + 0(1) as 
x — > oo. The quantity 

Q = b+l/b (18) 

is called the "background charge". The two-point function of the exponential of the 
boundary field behaves like 



^a(l>Biy)pa(j)Biy') 



Liouv 



d{a\z, zb) 

|y yi\2a{Q-a) 



(19) 



The explicit form of the function d was found in ref. fTO 



d{a\s) 



TTzT{b^)h- 



r(i - b 

X exp < - 



°o dt 

-oo t 



G{Q - 2a) 
G{2a - Q) 
sinh[(Q - 2a)t] cos^{st/b) (Q - 2a) 



InG(x) 



°° dt 
T 



:i-e-^*)(l-e-*/'') 



sinh(6t) sinh(6/t) 

(Q/2 - x)' 



+ 



t 



^ {Q/2 - x) 



t 



(20a) 
(20b) 
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where F denotes the Gamma function. The auxihary variable s depends on z and zb a.s 
follows 



cosh^fvrs) = — sinfTrfe^l 



(21) 



It is either real or pure imaginary. Writing 



s = u + iv 



u,v E R 



(22) 



one has namely 



^sin(7r62) > 1 
^sm(nb') < 1 



V = 0,u E {0, oo) 
u = 0,ve (0,1/2) 



(23a) 
(23b) 



The function d clearly fulfils the unitarity relation 



d{a\s)d{Q — a\s) = 1 



(24) 



It is straightforward to show from (jTU]) that the boundary exponential operators exhibit 
the following reflection property 



e'^'^s(?/) = d{a\s)e^^-''^'t'^{y) 



The associated boundary sinh-Gordon model is defined by the action 



(25) 



•^shci^, Zb) 



dV 



x>0 



— (V0)2 + 2zcosh(260) 

477 



/oo 
dy cosh(60B) (26) 
-oo 



The vacuum expectation value {e°''^'^)shG is conjectured [|] - [T^ to satisfy a reflection 
relation similar to (pSD 

(e'^^").hG = rf(ak)(e(«-"^'^^)shG (27) 
This relation, together with the obvious symmetry 



shG 



:28) 



determines the expectations up to a periodic function; the solution of interest is a "minimal 
solution" to the functional equations (^) and (pS)). 
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The boundary sine-Gordon model with the action ( |12b|) results as an analytical con- 
tinuation of (p^) via the substitutions b ^ ib, z ^ —z,zb —zb- The final result for 
{e^"''^^)sG = (e'"'^-^)z,^s reads [see eqs. (4.10) and (4.11) of ref. |]TD[, add a missing factor 
1/2 to ln^5(a)] 

'iizTil-b'^] 



where 

In (a) 

In gs{a) 



dt f 2sinh^(a6t) 
1 



2(1-6^) 



r(^2A 9o{a)gs{a) 
g(i-62)t/2 cosh(t/2) cosh(62t/2) - 1 



(29) 



sinhtsinh(62t) sinh((l — b'^)t) 
&t 2sinh^(a6t) sin2(st) 



a^e"* J>(30a) 



(30b) 



/o t sinhtsinh(6^t) sinh((l — 6^)t) 
where we have taken into account that s* = u — iv = ±s. This result holds under 
the conformal normalization of the bulk and boundary fields corresponding to the short 
distance asymptotics 

1 



I r — r^l^a^ 
1 



\y-¥ 



l\2a? 



as |r — t'I — > 
as \y — y'l ^ 



(31a) 
(31b) 



0, 0. This normalization is consistent with a well known leading short-distance be- 
haviour of the positive-negative pair correlation in the Coulomb gas. For the case of 
interest a = b, using the integral representation 



lnr(a;) 

one gets after some algebra 
9o{b) = 

9s{b) = 



°° dt 
— < 
t 



X-1 + 



1 



-{x-l)t _ I 

1 -e-* 



Re X > 



2v/^(l - 62)21(62) \2-2by \2-2b\ 
sinh 



'1-6" 



1-62 



sinh(7rs) 



(32) 

(33a) 
(33b) 



With the aid of (piD, one then finds 



47r3/2(l -62)2^ ^2 - 262^ \^2 - 262^ 



271 Zb 
f(62 



n 1/(1-62) 



2 cosh(7rs) 



1-"') sinh (^) 
sinh(7rs) 



(34) 
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According to (|21|), in the limit z ^ the real parameter s diverges and (0) yields 



1 



in full agreement with formula (24) in ref. M. 



2ttzb 
f(62) 



l/{l-f,2) 



(35) 



4 Surface tension 

We first summarize the bulk thermodynamics of the 2D Coulomb gas The pressure 
is expressible in terms of the soliton mass M as follows 



^■^4 V 2 



(36) 



where q = /3/{A — /S). The normalization of the bulk field ( |31a| ) fixes the relationship 
between the fugacity z and M in the form 



r(g/(g + i)) 



M 



v^r((g+l)/2) 



1 2/(g+l) 



7rr(l/(g+l)) [ 2r(g/2) 
The total particle number density n, generated via 

d 



(37) 



is related to M as follows 



n 



n = z—pp 
oz 



lM2(l + g)tan(^ 



(3^ 



(39) 



Let us denote the surface tensions of the metal-electrolyte and of the dielectric- 
electrolyte boundaries by 



Jmet = hm 'j{z,Zb) 

ZB^OO 



(40a) 
(40b) 



respectively. In paper |^ [formula (33b) with zero bulk potential ^ = 0], we have found 

M 



met 



4cos(g7r/2) 



sm 1 Y I - cos [ — ] + 1 



(41) 



11 



With regard to the definitions (|40|) and eq. (0), it holds 

d 



Pi, 



diel 



met 



B- 



-h{z,ZB) 



dzB 



Z=0,Zg 



(42) 



Inserting the previously obtained formulae (|3^ ) and (^) into ( ^21) and substituting b 
V^/2 in accordance with ( |llc| ), one arrives at 



Pldiel = Pimet + 



7r3/2(4-/5) 



2tt 



r(/3/4) 



4/(4-/3) 



r 



'2-/3' 



r 



2/5, 



[/i + h) (43) 



where 



/i 



2n 



1 2/(4-/3) 



4sin(7r/3/4) 



dw < sinh 



2tt 



I Atxu > 

-, 2/(4-/3) 



— sinh(7rn) [2 cosh(7rn)] 



/3/(4-/3) 



(44a) 



4sin(7r/5/4) 

1/2 r / 

dt" < sin 



At^v 



sin(7rf ) [2 cos(7rt')] 



/3/(4-/3) 



(44b) 



Here, by using eq. (pT]) we have transformed the integration over Zb into integrations over 
the M, f -components of the auxiliary variable s (^2]); the split of the integral onto the /i 
and I2 ones is due to the existence of the two regimes (|23a|) and ( p3b|) . Algebra gives 



h 
h 



M-P) 



2(4-/5) 



4sin(7r/?/4) 

z 



2/(4-/3) 



1^2/3/(4-/3) _ 



-I 2/(4-/3) 



4sin(7r/5/4) 



^ _ 2/3/(4-/3) 



cos 



^ 27r ^ 



(45a) 
(45b) 



Inserting Ji and /2 into (^3]), then using (|37|), applying simple operations with trigono- 
metric functions and formula r(a;)r(l — x) = n / sin(7rx) and considering (|41|), one finally 
gets 



Pi, 



M 



diel 



sm 



cos 



4cos(g7r/2) 

where q = P/{A — jd) and M is related to the fugacity z by 



1 



(46) 
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When M is expressed in terms of the particle density n by using eq. (|39D, (^61) takes 
the form 



diel 



n(4 - (3) 



2sin(7r/3/(4-/3)) 



1/2 



sm I — I + cos 



2(4-/3), 



2(4-/3), 



- 1 



(47) 



P^diei has the small /3-expansion 



1/2 



(4J 



The singular behaviour oi n as (3 ^ 2 {z fixed) can be deduced from formulae (^) and 



n 



2^ 



(49) 



On the other hand, the surface tension, expressed in terms of the fugacity z, remains 
finite at the bulk collapse point (3 = 2, 



irz 



Idiel = 13 = 2 



(50) 



The function jS'-^diei increases monotonously up to /3 = 3 where the surface tension diverges 

-r(i/4)-^ 



TT 



r(3/4) 



3-/3' 



/? ^ 3- 



(51) 



The singularity prevents from going beyond this point. 



5 Analytic checks 

The surface tension 7 (= jdiei) is the boundary part per unit length of the grand potential 
Q. The total number of particles is given by = A^_|_ + = —(3zdQ/dz. The boundary 
part of this relation is 

d'-f f°° 

— (3z— = dx [n{x) — n] (52) 
oz Jo 

or, since z oc n^~^/^, 

- (3n^ = 1^1 - ^ j dx [n{x) - n] (53) 
These formulae will be used for computing the surface tension from the density profile. 
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5.1 High-temperature expansion 



As a check of the exact expression (^T]) for the surface tension, its high-temperature 
expansion in (3, equation (^81) , can be compared to a direct evaluation of the first two 
terms of this expansion by using a renormahzed Mayer expansion, in close analogy with 
section 5 of paper [^]). Since the derivative of the interaction potential f (r, r^) +v{r* , r') 
with respect to x vanishes when point r is on the interface, the same boundary condition 
holds for the renormahzed bond K{r,'r'), 



After some algebra, one finds 



/ dx [n{x) — n] 
Jo 



1 1/2 



16 



_1 - f (iH- ■ 



Consideration of ( p5D in (|53D gives the final result 



TT 



in agreement with 



(54) 



(55) 



(56) 



5.2 Free fermion /3 = 2 case 

Using a Thirring-like representation of the Coulomb gas in presence of a dielectric wall at 
special inverse temperature (3 = 2, the exact result for the density profile was obtained in 



the form |ri| 



n±{x) —n± = ——KQ{2mx) 



(57) 



where m = 2'kz is the rescaled fugacity. Thus, 

/•oo 

/ dx[n(x) — 
Jo 



TTZ 



n\ 



(5^ 



Inserting this integral into equation (j5^), one immediately arrives at the exact formula 
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6 Concluding remarks 

A two-dimensional model for the interface, the two-component plasma bounded by a rec- 
tilinear ideal dielectric wall, was considered. The surface tension as a function of fugacity 
is finite at any inverse temperature /5 < 3. In general, the surface properties are supposed 
to be governed by the particle-image interaction. A particle at distance x from the di- 
electric wall interacts with its own image of the same sign through a repulsive potential 
— (1/2) ln(2a;), and the corresponding Boltzmann factor (2a;)^/^ is always integrable at 
small X, and gives rise to a vanishing particle density at the interface. In spite of this, our 
exact solution for the surface tension predicts a surface collapse at /3 = 3, identified with 
the divergence of the surface tension with a characteristic type of short-distance collapse 
singularity, see formula (^l]). Such a phenomenon might be explained by recalling an 



observation of Hansen and Viot [^] concerning the short-distance behavior of bulk pair 
distribution functions for two charges of the same, let us say plus, sign. The expected 
dependence ~ C'^+r^ as r — > 0, where the pref actor C++ is related to a difference of 
free energies, is changed to (?++ ~ C+^r^^^ at /9 > 1 as a consequence of the divergence 
of C++ at point (3 = 1. The weaking of (?++ at short distance is caused by a pair forma- 
tion of oppositely charged particles at low temperatures: the neutrality of a pair allows 
a third particle to approach very close to the pair. For (3 > 2, the strong clustering of 
positive-negative particles causes an effective short-distance attraction between particles 
of the same sign. Based on these plausible arguments it is tempting to conjecture that, 
at /3 = 3, the surface tension diverges due to a paradoxical short-distance collapse of a 
particle with its own image of the same sign. 
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